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On the spectrum of the pulsed gamma-ray emission from 
lOMeV to 400GeV of the Crab pulsar. 
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ABSTRACT 

In the present paper a self-consistent theory, interpreting the VERITAS and the 
MAGIC observations of the very high energy pulsed emission from the Crab pulsar is 
considered. The photon spectrum between lOMeV and 400GeV can be described by 
two power-law functions with the spectral indexes equal to 2 and 3.8. The source of 
the pulsed emission above lOMeV is assumed to be the synchrotron radiation, which is 
generated near the light cylinder during the quasi-linear stage of the cyclotron insta- 
bility. The emitting particles are the primary beam electrons with the Lorentz factors 
up to 10 9 . Such high energies by beam particles is supposed to be reached due to 
Landau damping of the Langmuir waves in the light cylinder region. This mechanism 
provides simultaneous generation of low (radio) and high energy (10MeV-400GeV) 
emission on the light cylinder scales, in one location of the pulsar magnetosphere. 
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1 INTRODUCTION 

The recent observations of the Crab pulsar in the very 
high energy (VHE) domain with the VERITAS array of 
atmospheric Che renkov telescope s revealed pulsed 7-rays 
above lOOGeV (|Aliu et all l201lft . which was later con- 
firmed by measurments of the MAGIC Cherenko v telescope 
' Aleksic et all l201ll . 12012ft . Prior to the work lAliu et all 



2011ft the highest energy of the pulsed em ission from the 



Crab pulsar was 25GeV (|Aliu et all 120081 '). The detection 
of such a high energy pulsed 7-rays cannot be explained 
on the basis of current pulsar emission models. It is gener- 
ally assumed that the VHE emission is produced either by 
the Inverse Compton scattering or by the curvature radi- 
ation. By analyzing the aforementi oned emission processes 
|Machabeli fc Osmanovll2009l . [2?5lol '). we have found that for 
Crab pulsar's magnetospheric parameters even very ener- 
getic electrons are unable to produce the photon energies 
of the order of 25GeV. Studying the curvature radiation, it 
was shown that the curvature drift instability efficiently rec- 
tifies the magnetic field lines, leading to a negligible role of 
the curvature emission process in the obser ved VHE domain 
Osmanov et all 120091 ). In previous work Chk heidze et al.l 



2011ft we have explained the origin and the measured spec- 



trum of the Crab pulsar in the high energy (HE) domain 
(0.01 — 25)GeV, relying on the pulsar emission model first 
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developed bv lMachabeli fc Usovl (|l979l h According to these 
works, in the electron-positron plasma of a pulsar magne- 
tosphere the low frequency cyclotron modes, on the quasi- 
linear evolution stage create conditions for generation of the 
HE synchrotron radiation. 

It is well known that the distribution function of rela- 
tivistic particles is one dimensional at the pulsar surface, be- 
cause any transverse momenta (p±) of relativistic electrons 
are lost in a very short time (^ 10 _20 s) via synchrotron 
emission in very strong magnetic fields. But plasma with an 
anisotropic one-dimensional distribution function is unsta- 
ble which inevitably leads to the wave excitation process. 
The main mechanism of the wave generation in plasmas of 
the pulsar magnetosphere is the cyclotron instability, which 
develops at the light cylinder length-scales (a hypotheti- 
cal zone, where the linear velocity of rigid rotation exactly 
equals the speed of light). During the quasi-linear stage of 
the instability a diffusion of particles arises as along, also 
across the magnetic field lines. Therefore, plasma particles 
acquire transverse momenta and, as a result, start to radiate 
in the synchrotron regime. 

In IChkheidze et alj (|201ll h it was shown that near the 
light cylinder the radio waves are generated, provoking the 
re-creation of the pitch angles and the subsequent syn- 
chrotron radiation in the HE domain. Thus, in the frame- 
work of this model generation of low and high frequency 
waves is a simultaneous process and it takes place in one 
location of the pulsar magnetosphere. This explains the ob- 
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served coincidence of the HE emission pulses with the ra- 
dio signals. The recent VERITAS and MAGIC observations 
have shown that the aforementioned coinci dence takes places 
in th e VHE domain (> l OOGeV) as well jAliu et al.ll2008l . 
l201ll ; lAleksic et alJlioITh . According to lAliu et alj (|201lh 
detection of pulsed 7-ray emission of the order of lOOGeV 
requires that the emission should be produced far out in the 
magnetosphere. Thus, we suppose that the pulsed high and 
the very high energy radiation of the Crab pulsar is gener- 
ated through the synchrotron mechanism at the light cylin- 
der length-scales, switched on due to the quasi-linear diffu- 
sion. The resonant particles are the primary beam electrons 
with the Lorentz-factor 76 ~ 10 8-9 , giving the synchrotron 
emission in the ( 0.01 — 400)GeV e nergy domain. 

According to lAliu et ail (|2008h a joint fit to the EGRET 
(lOMeV to lOGeV) and MAGIC (> 25GeV) data predicted 
a power-law spectrum with a generalized exponential shape 
for the cutoff, described as F e oc e~ a exp(— (e/eo)' 9 ), where 
q = 2.022 ± 0.014. We provided a theoretical confirmation 
of the measured spectrum, which yielded j3 = 1.6 and the 
cutoff energy e = 23GeV (|Chkheidze et al.l 120111 ). Recent 
VERITAS observations (100 - 400) GeV combined with the 
Fermi-LAT data (0.1 — 10) GeV favor a broken power law as 
a parametrization of the spectral shape. The good fit results 
are also obtained if one uses a log-parabola function, but it 
fails to describe the spectrum below 500MeV. Although, the 
Fermi-LAT and Magic data below 60 GeV can be equally 
well parameterized by broken power law and exponential 
cutoff. In the energy range between lOOGeV and 400GeV 
measured by VERITAS and MAGIC, the spectrum is well 
described by a simple power law with the spec tral index 
equal to 3.8 jAliu et al.ll201ll ; lAleksic et alj|2012h . 

The paper is organized as follows. In Sect. 2 we describe 
the emission model, in Sect. 3 we derive the theoretical syn- 
chrotron spectrum for the high and the very high energy 
7-ray emission of the Crab pulsar and in Sect. 4 we discuss 
our results. 



electron's charge and the rest mass, c is the speed of light, 
k x is the wave vector's component along the drift and B 
is the magnetic field induction. During the wave generation 
process, one also has a sim ultaneous feedback o f these waves 
on the resonant electrons (|Vedenov et al.|[l96ll ). This mech- 
anism is described by the quasi- linear diffusion (QLD), lead- 
ing to a diffusion of particles as along as across the magnetic 
field lines. Therefore, resonant particles acquire transverse 
momenta (pitch angles) and, as a result, start to radiate 
through the synchrotron mechanism. 

The wave excitation leads to redistribution process of 
the resonant particles via the QLD. The kinetic equation 
for the distrib ution function of the resonant electrons can 
be written as |Chkheidze et al1l201lh : 



df° 
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P± op± [ dp± 
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are the transversal and longitudinal components of the syn- 
chrotron radiation reaction force, where a a — 2e 2 oj%/3c 2 
and D± 1 is the transverse diffusion coefficient which is de- 



fined as follows ijChkheidze et al.ll201ll ) 



8mccj^7| 



\E, 



(4) 



and its value can be estimated from the expression \Ek\ ~ 
mc 2 nt,7f,c/2tj c , where uj c is the frequency of the cyclotron 
waves. From Eq. (1) it follows that 



LOB 



(5) 



2 EMISSION MODEL 

Any well known theory of pulsar emission suggests that, 
the observed radiation is generated due to processes tak- 
ing place in the electron-positron plasma. It is generally as- 
sumed that the pulsar magnetosphere is filled by dense rela- 
tivistic electron-positron plasma with an anisotropi c one- 
dimen sional distribution function (see Fig. 1 from lArons! 
(1981)) and consists of the following components: a bulk of 
plasma with an average Lorentz-factor 7 ~ y p , a tail on the 
distribution function with 7 ~ ft, and the primary beam 
with 7 ~ 7t. The distribution function is one-dimensional 
and anisotropic and plasma becomes unstable, which might 
cause a wave excitation in the pulsar magnetosphere. The 
main mechanism of wave generation in plasmas of the pul- 
sar magnetosphere is the cyclotron inst ability. The cyclotro n 
resonance condition can be written as (|Kazbegi et al.|[l992h : 



k„V„ 



k x u x + — = 0, 

It 



(1) 



where u x = cV^ r / puiB is the drift velocity of the particles 
due to curvature of the field lines with the curvature radius 
p, lob = eB/mc is the cyclotron frequency, e and m are the 



where 8 = lo 2 / '(4^%^) , co p = 4Tvn p e 2 /m is the plasma 
frequency and n p is the plasma density. 

The transversal QLD increases the pitch-angle, whereas 
force Fl resists this process, leading to a stationary state 
(df/dt — 0). The pitch-angles acquired by resonant electrons 
during the process of the QLD satisfies ^ — p±/p\\ -C L 
Thus, one can assume that d/dp± >> d/dp\\. In this 
case the solution of Eq.(2) gives the distribution function 
of the resonant partic les by their perpendicular momenta 
l|Chkheidze et alJl201lD 



flj>±) = Cexp 



-dp j 
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where 



n 1 / 2 /3m 9 c n 7 6 5xl/4 
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And for the mean value of the pitch angle we find ipo ~ 
p± /p\\ — 10 -6 . Synchrotron emission is generated as the 
result of appearance of pitch angles. 

The synchrotron emission flux of the set of electrons in 
the framework of the presen t emission scenario is written as 
(see lChkheidze et all (|201ll )) 
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(8) 



Here / (pj| ) is the longitudinal distribution function of elec- 
trons, e m ~ 5 • 10 _18 -B?/'o7 2 GeV is the photon energy of the 
maximum of synchrotron spectrum of a single electron and 
K$/ 3 (z) is the Macdonald function. After substituting the 
mean value of the pitch-angle in the above expression for 
e m , we get 



5 ■ 10" 
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Accordingly, the beam electrons should have 75 ~ (6 ■ 
10 8 — 10 9 ) to radiate the photons in the energy domain 
~ (10 — 100) GeV energy. The gap models provide the 
Lorentz factors up to 10 7 , which is not enough to explain 
the detected pulsed emission. Consequently, the additional 
particle acceleration mechanism should be invoked to accel- 
erate the fastest electrons to even higher energies. 

According to lAliu et all l|201ll ) the observations of the 
Crab pulsar indicates that the HE pulsed emission should 
be produced far out in the magnetosphere. Therefore, we 
assume that the Langmuir waves (L) generated via the 
two-stream instability at the light cylinder length-scales 
undergo Landau damping on the fastest beam electrons, 
which results in their effective acceleration. Excitation of 
L waves through the two-stream instability in the rela- 
tivistic electron-p ositron plasma is considered in a series 
of wor k s (see e.g.. lLominadze fe Mikhailovskiil (Il979l ) ; lUsovl 



(199S 



Ursov fe Usov (1198. 



(1 19871); lAsseo fc Melikidzd 

iGogoberidze et al.l (|2008l V), where it is assumed that the in- 
stability develops due to overlapping of the fast and slow 
pair plasma particles at distances r ~ 10 8 cm from the star 
surface. For typical parameters of pair plasma in the pulsar 
magnetospeher s the growth rate o f the instability is quite 
sufficient (see IGogoberidze et ah! |2008t )). that inevitably 
provides existence of L waves in the vicinity of the light 
cylinder zone. The phase velocity of the excited L waves is 
asym ptotically close to the speed of light l|Gogoberidze et ah! 
2008). Therefore, these waves can be only damped on the 
fastest electrons of the primary beam, which velocity equals 
the phase- velocity of the waves and the distribution function 
satisfies df bo /dp\\ < 0. Taking into account the equipartition 
of energy among the plasma components n P 7 P = n b ^ b /2, 
one can estimate the energy density of L waves as the 
half of the energy density of the primary beam particles. 
Through the Landau damping process the energy of the 
waves is transferred to the small fraction of the beam elec- 
trons with the highest Lorentz factors. If we assign the 
density of these electrons as n* and the Lorentz factors as 
7*, and equate the energy densities of particles and the L 
waves, we will get 2n* /rib ~ lb/"/*- The total density of 
the beam electrons is equal t o the Goldreich- Julian de nsity 
rib = B/Pce « 2 • 10 7 cm" 3 jGoldreich fe Julianlfl969l ) and 
also if we take into account that ji, ~ 10 7 and 7* ~ 10 9 , 
we find n* ~ 10 5 cm" 3 . As we see, if the wave energy is 
transferred to the fastest beam electrons, which number is 
two orders of magnitude smaller than the total number of 
the primary beam particles, they will gain the Lorentz fac- 
tors up to 10 9 . The observational fact, that the emission flux 



above 25GeV decreases should be caused by reduced number 
of emitting particles with the highest Lorentz factors. 

It should be mentioned that during the Landau damp- 
ing process the beam distribution function will be alongated 
and will form a high energy 'tail'on the distribution function. 
The final shape of the distribution function after the quasi- 
linear relaxation is the plateau, and the stationary state is 
reached. But in our case this might not be achieved as in the 
same region where the L waves are damped the cyclotron in- 
stability is developed, which involves the beam electrons into 
the cyclotron resonance process. This complicated process 
causing the redistribution of the resonant particles needs a 
more detailed investigation, which we plan to perform in 
our future work. In the present paper, we estimate the final 
shape of the distribution function which provides explana- 
tion of the measured spectrum. 

The beam particles lose their energy through the syn- 
chrotron radiation, which sets the upper limit on the Lorentz 
facto rs that can be achiev ed during their acceleration pro- 
cess (|de Jager et al. II 19961 ). The maximum achievable value 
of 7 can be estimated by equating the synchrotron radiative 
looses to the power of the emitting particles gained through 
the acceleration process, which in our case can be written 
as: 



2e*BW 2 p 

— = mc ~/Fld- 

3 ra A c 6 

Here 



5/2 ' 



(10) 
(11) 



is the Landau d amping rate, where uj b — <J ' A-Ke 2 n b /m 
l|Volokitin et al. I Il987l ). The estimations show, that the 
Lorentz factor of the beam electrons that can be reached 
through the acceleration process 7 < 10 19 . Consequently, 
the limit on the energy of synchrotron photons in our case 
should be e ma x ~ 10 54 eV (see Eq. (9)), which inevitably al- 
lows generation of the observed ~ lOOGeV photons from the 
Crab pulsar through the synchrotron emission mechanism. 



3 SPECTRUM OF THE SYNCHROTRON 
RADIATION 

To obtain the synchrotron emission spectrum in our case, 
we need to solve the integral (8). For this reason, first let us 
find the parallel distribution function of the beam electrons 
L during the QLD process of the cyclotron instability. By 
multiplying both sides of Eq. (2) on p± , integrating it over 
p± and taking into account that the distribution function 
vanishes at the boundaries of integration, Eq. (2) reduces to 



dt 



d 

dpu 



o/||) 



(12) 



For jip -C 10 10 , a magnetic field inhomogeneity does 
not affect the process of wave excitation. The equation that 
describes the cyclotron n oise level, in this case, has the form 
l|Lominadze et al~lll983h 



9\E k \' 2 . |2 

where 



(13) 
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is the growth rate of the instability. Here fey can be found 
from the resonance condition (1) 



cSj r 



Combining Eqs. (10) and (11) one finds 
= 0, 
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Consequently, one can write 
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d 

"apu 



1/2 



const. 



(15) 



(16) 



(17) 



(18) 



Taking into account that for the initial moment (the moment 
when the cyclotron instability arises) the energy density of 
cyclotron waves equals zero, the corresponding expression 
writes as 



d 

dp\\ 



1/2 



■f\\i 



(19) 



Let us assume that \E k \ oc i _m (as there is no direct way 
to calculate the dependence |£Ife(-y)|, we can only make an 
assumption and check its plausibility by fitting the theoret- 
ical emission spectrum with the observed one), in this case 
for the parallel distribution function we will have 



h K 4, +7~ 



(20) 



The initial distribution f» of the beam particles in this case 
is the redistributed one after the Landau damping process. 
The final shape of the distribution function of resonant par- 
ticles via the Landau damping when the stationary state is 
reached is the plateau. But in our case this might not be 
reached as in the same region develops the cyclotron insta- 
bility. Thus we consider /j, oc 7~™, where n is not close 
to zero. Consequently, the parallel distribution of the beam 
electrons is proportional to two power-law function with the 
indexes m + 3/2 and n. 

The e ffective value of the pit ch angle depends on \Ek\ 2 
as follows l|Chkheidze et alj|201lh 
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(21) 



Using expressions (8), (18) and (19), and replacing the inte- 
gration variable p« by x — e/e m , we will get the synchrotron 
emission spectrum 
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where 
Gi(y)- 



K 5 / 3 (z)dz 



dx 



G 2 (y) 



K 5/3 (z)dz 



dx. 



(23) 



10 y in which case, we have (e/e m )max -C 1 
mjrmn 3> 1. Under such conditions the functions 



The energy of the beam particles vary in a broad range 
76 ~ 10 6 1 " 9 
and (e/e. 

Gi{y) ~ G-2(y) ~ G(0). Consequently, we can assume that 
the synchrotron emission spectrum (Eq. (20)) is propor- 

_ 2ml4n-l _ 6m + 5 

tional to two power-law functions e 5-2™ anc [ € 5-2™ _ 
According to VERITAS observations the spectrum mea- 
sured in the energy domain (100 — 400) GeV is wel l described 
by po wer-law with the spectral index equal to 3.8 (jAliu et al.l 
l201lh . When m « 1, we have -(6m + 5)/(5 - 2m) = -3.8. 
At the same time, when m = 1 and n — 1.2, we find 
— (2m + 4n — l)/(5 — 2m) = —2 that is in a good agreement 
with the observations in the lOMeV - 25GeV energy domain, 
which shows the power-law spectrum F(e) oc e - 2 022 ± - 014 
ijAliu et alJl2008h . 



4 DISCUSSION 

The interesting observational feature of the Crab pul- 
sar of the coincidence of pulse-phases from different fre- 
quency bands, ranging from radio to VHE gam ma-rays 
^Manchester fc Tavlorlfl980l ; lAliu etHI 120081 . [ioilh implies 
that generation of these waves occur in one location of 
the pulsar magnetosphere. This consideration automatically 
excludes the generally accepted HE emission mechanisms, 
the Inverse Compton (IC) scatter ing and the curvature 
radia t ion, w hich are not localized l)Machabeli fc Osmanovl 
120091, |201Ch. These particu lar issues have been considered 
by Machabeli fc Osmanovl (|2010h . Studying the curvature 
radiation it was found that the curvature drift instability is 
efficient enough to rectify the magnetic field lines (curvature 
tends to zero) in the region of the generation of high and the 
VHE emission, making the curvature emission process neg- 
ligible. At the same time by analyzing the IC scattering, it 
was found that for Crab pulsar's magnetospheric parame- 
ters even very energetic electrons are unable to produce the 

observ ed HE photons. 

In lLvutikov et all l|201lf ) it is argued that the main gen- 
eration mechanism of the VHE emission of the Crab pulsar is 
the IC scattering of the soft UV photons by the secondary 
plasma particles. Particularly, it is assumed that the sec- 
ondary plasma, which is produced due to cascades in the 
outer gaps of the magnetosphere is responsible for the soft 
UV emission via the synchrotron mechanism. This radiation 
plays the target field role for the IC scattering process. As 
a result, the VHE 7-ray emission is produced extending to 
hundreds of GeV. Let us consider the equation that gives 
the frequency of the phot on after the IC scattering (e.g. 
iRvbicki fc Light manl l|l979h ): 



(1 - /3cos6) 



(24) 



1 - /3cos0' + (1 - cosd")hu}/jmc 2 ' 

where ui is the frequency before scattering, /3 = v/c, 9 = 
(PK), 0' = (PK 7 ) and 6" = (KK) (by P we denoted the 
momentum of relativistic electrons before scattering, K and 
K' denotes the three momentum of photon before and af- 
ter scattering, respectively). The emission maximum comes 
along the magnetic field lines. The pitch angles of the rela- 
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tivistic electrons moving along the magnetic field lines are 
very small (ip ~ 10 -6 see Eq.(7)), accordingly as the UV 
emission is generated through the synchrotron regime, the 
angle 9 should also be very small. Since we observe the well- 
localized pulses of the VHE emission, the angle 9' <C 1. Tak- 
ing into account the observational fact of the phase coinci- 
dence of signals from different frequency bands, one should 
assume that 6 ~ 9'. At the same time the coincidence of 
pulse signals of UV and the VHE emission gives 9" = 0. 
Consequently, the IC scattering of the soft UV photons by 
the secondary plasma electrons in case of the Crab pulsar 
should only cause the redistribution of the electrons, but can 
not provide the increase in the emission frequency, especially 
up to the VHE band. For significantly increasing the photon 
energy by the IC scattering processes without violating the 
condition of the pulse-phase coincidence, the angle 9 must 
be large enough, which can not be provided by our model. 
Apparently, the IC scattering should play the main role in 
the generation of the high energy emission for pulsars with 
the large pulse profiles in soft energy domains. 

The emission model proposed in the present paper im- 
plicitly explains the observed pulse-phase coincidence of low 
(radio) and high frequency (10MeV-400GeV) waves, as their 
generation is a simultaneous process and it takes place in 
the sa me place of the pulsar magnetosphere. In previou s 
works (|Machabeli fc Osmanovll20ld ; IChkheidze et al.ll201ll ') 
we applied this model to explain the HE (0.01-25GeV) 
pulsed emission of the Crab pulsar obse rved by the MAGIC 
Cherenkov Telescope l|Aliu et al.l 120081 ). It was found that 
on the light cylinder length-scales the cyclotron instability 
is arisen, which on the quasi-linear stage of the evolution 
causes re-creation of the pitch angles and as the result the 
synchrotron radiation mechanism is switched on. We assume 
that the source of the high and the VHE pulsed emission of 
the Crab pulsar is the synchrotron radiation of the ultrarel- 
ativistic primary beam electrons. To explain the observed 
high frequency gamma-rays by synchrotron mechanism the 
Lorentz factors of the emitting particles should be of the or- 
der of 10 s - 10 9 (see Eq. (9)). The highest Lorentz factor for 
the typical pulsar is ~ 10 7 . Thus, we assume that such an 
effective particle acceleration is caused by existence of the 
Langmuir waves with the phase velocities v v < c close to the 
velocities of the fastest beam electrons in the region close to 
the light cylinder. Consequently, the L (electrostatic) waves 
are efficiently damped on the most energetic primary beam 
electrons. The Landau damping causes the growth of a HE 
tail on the distribution function of the resonant electrons 
and inevitably throws the most energetic particles to higher 
Lorentz factors (up to 7 ~ 10 9 ). At the same time the beam 
electrons acquire the pitch angles due to the cyclotron in- 
teraction with the transverse waves, which causes the syn- 
chrotron radiation processes giving the observed high and 
the VHE emission up to 400GeV. The distribution func- 
tion tends to form the plateau (due to Landau damping), 
though the number of processes impede this. The reaction 
force of the synchrotron emission, scattering of L waves and 
the Compton scattering of photons on the beam particles 
also take place in process of formation of the distribution 
function of beam electrons. As the result, it is unlikely for 
the distribution to reach the shape of plateau and thus, we 
represent it as fu oc 7"". 

The calculation of the synchrotron emission spectrum 



by taking into account the processes described above (see 
Eq. (20)) and matching it with the observations shows that 
the emission spectrum in the 10MeV-25GeV energy domain 
depends on the power-law index n of distribution function of 
the beam particles. At the same time the emission spectrum 
in (100 — 400) GeV energy domain does not depend on n but 
only depends on m (see Eq. (18)). When m = 1 and n = 1.2 
the emission spectrum well matches the measured one in 
both high (0.01 - 25GeV) and the ve ry high (100-400GeV) 
energy domains. In previous paper l|Chkheidze et al1l201ll ) 
explaining the HE (10MeV-25GeV) spectrum we obtained 
the power-law function with the exponential cutoff F e oc 
€~ 2 exp[— (e/23) 1 ' 6 ], as jb ~ 10 s were the highest considered 
Lorentz factors of the emitting electrons. We assume that 
detection of the exponential cutoff at higher (> 400GeU) 
photon energies can not be excluded, the exact location of 
the cutoff energy is defined by the highest energy of the 
beam electrons that can be reached through the Landau 
damping before the particles reach the light cylinder. This 
particular problem needs more detailed investigation, which 
is the topic of our future work. 
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